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Abstract 

<N 

An exponential-type approximation of the first order Marcum Q-function is presented, which is robust to changes 
in its first argument and can easily be integrated with respect to the second argument. Such characteristics are 
particularly useful in network connectivity analysis. The proposed approximation is exact in the limit of small first 
| argument of the Marcum Q-function, in which case the optimal parameters can be obtained analytically. For larger 

values of the first argument, an optimization problem is solved, and the parameters can be accurately represented 
| using regression analysis. Numerical results indicate that the proposed methods result in approximations very close 

to the actual Marcum Q-function for small and moderate values of the first argument. We demonstrate the accuracy 
of the approximation by using it to analyze the connectivity properties of random ad hoc networks operating in a 
Rician fading environment. 

I. Introduction 
^ ' The Marcum Q-function, defined as the integral CQ 



Qi(a,b) = I xexp ( ) Io(ax)dx (1) 



> 

\Q . for a, b > where Iq(x) is the modified Bessel function of the first kind, is a fundamental function that arises 

O I in the performance evaluation of a wide class of communication systems CQ-lH]]. From a mathematical point 

* • ■ of view, this function represents the complementary cumulative distribution function (CCDF) of the power of a 

^ . Rician distribution. The integral representation of the function given by © cannot be manipulated easily to provide 

£S) [ simple expressions for the performance of communication systems, especially when the function Qi(a, b) must be 

" integrated with respect to one of its arguments (T|. To solve this issue, numerous works have proposed alternative 

. representations of Qi(a, b) to facilitate analysis (see, e.g., O, H and references therein). Exponential-type bounds, 

! provided they are tight, have been particularly attractive, especially when evaluating the bit error rate at high 



signal to noise ratio (SNR) (3J, O. In other situations, approximations may be more suitable than bounds [6] 



0. However, such approximations may have complicated mathematical structures and/or be inaccurate in certain 
domains of their arguments. 

In this paper, a simple exponential approximation of the first order Marcum Q-function is presented that yields 
small approximation error over a large domain in its two arguments. The approximation is designed such that it 
can be used in situations where Q±(a, b) must be integrated over its second argument. In what follows, a heuristic 
approach is first employed to find the right form of the approximation, which is parameterized by two functions 
of a. An analytical framework for determining the correct parameterization is then explored, which is shown to 
be accurate for < a -C 1. For a 3> 1, the optimal parameterization is calculated numerically. Although the 
proposed approximation is useful in its own right, it is particularly helpful in situations where the Q-function 
must be integrated. To illustrate this fact, we present an example application of our results whereby the proposed 
approximation is used to analyze the connectivity probability of a random ad hoc network operating in Rician 
fading channels. 
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The structure of this paper is as follows. In Section [TTJ the proposed approximation and means of deriving the 
optimal a-dependent parameters are presented. An example application of the approximation is given in Section ITTTl 
while the accuracy of the approximation is presented in Section [IVJ Finally, some concluding remarks are given in 
Section |VJ 

II. Approximation of Qi(a, b) 

By plotting Q\(a,b) as a function of b for various values of a, it can be readily observed that Q\(a,b) decays 
exponentially with b, where the value of a roughly defines the shift of Q\ along the 6-axis. Consequently, we 
propose to approximate Q\{a,b) by the function 

Qi(a,6) = exp (-e^ a V (a) ) (2) 

where v{a) and //(a) are nonnegative parameters dependent upon a. The key is to choose these parameters such 
that the accuracy of the approximation is high. 

As previously discussed, we are concerned with obtaining an approximation that is useful over the range of the 
argument b for some fixed a. Thus, we define the approximation error as the function 

POO 

£(a)= / (QiM)-Qi(a,6)) 2 d&. (3) 
J o 

Furthermore, we define fi(a) and v{a) to be polynomials of order m, where a larger value of m yields a better 
approximation. Thus, we have 

fi(a) = no + fi\a + ^ 2 a 2 H h [i m a m 

v(a) = u + v\a + u 2 a 2 H 1- v m a m 

in which case the approximation becomes 

Qx(a,b) = exp f-eZn= (Hninb+v n )a"\ (4) 



The goal is now to choose the coefficients {fiQ, . . . , ju m , vq, . . . , u m }, independent of b, such that £(a) is minimized. 
Depending on the value of the argument a, this can be done analytically or numerically. 

A. Analytical Approach for Small Arguments 

First, consider the case where < a <C 1. It is logical to expand Qi(a, b) and Qi(a, b) about a = and equate 
the coefficients term by term. Of course, if the expansions for Q and Q converge and the corresponding coefficients 
match to arbitrary order, then Q = Q. Since this is clearly not the case, it is advisable to equate coefficients 
recursively, from lowest to highest order. 

For example, let m = 4. To leading order, we have Qi(0,6) = exp(— b 2 /2) and Qi(0,b) = exp(— e Ua b^°). It 
follows that we should choose /j,q = 2 and vq = — In 2 since this ensures the approximation is exact at a = 0. 
Next, we can equate the first order terms to obtain the equatiorQ [i\ In b + v\ = 0. But we see from (|4]) that this 
formula implies there is no 0(a) term in the second exponent of Q. Thus, we may take p,\ = v\ = to maintain 
independence of b. This process can be continued in a straightforward manner. However, when we equate the fourth 
order terms, we obtain the equation /^lnfe + z^ = 6 2 /32, and thus either ^4 or ^4 is dependent upon b, a condition 
that is not allowed by our definition of the polynomials \i and v. Instead, we can optimize £(0) over p,^ and ^4 by 
differentiating with respect to each variable, setting the results to zero, and solving for /X4 and ^4. This yields the 
optimal fourth order polynomials 

9 

u(a) = 2 H — ; — -a 

PV ' 8(9^ 2 -80) 

a? 45vr 2 + 72 In 2 + 36C - 496 4 

via) = — In 2 1 ; — Ti ; a 

v ' 2 64(9vr 2 - 80) 

which are independent of b, and thus satisfy the conditions of our approximatiorH By substituting these expressions 
for n(a) and v{a) into Q and evaluating the integral in © for small a, it is apparent that £(a) 7.5 x 10~ 5 a 8 . 
Thus, the fourth order result is very accurate for a <C 1, an observation that is corroborated by Fig. Q] 

'The details of the calculations are straightforward but lengthy, and are thus omitted here for brevity. 
2 The symbol C denotes the Euler-Mascheroni constant, where C ~ 0.5772. 
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Fig. 1. Approximation error derived from fourth order polynomial representation of fi(a) and v(a) for small a. 



TABLE I 

Suitable v and /i for different a 



a 


V 


i 1 


1.0000 


-1.1739 


2.0921 


2.0000 


-2.5492 


2.7094 


3.0000 


-4.6291 


3.6888 


4.0000 


-7.1668 


4.7779 


5.0000 


-10.0339 


5.9074 


6.0000 


-13.2014 


7.0794 



B. Numerical Approach for General Arguments 

While a closed-form expression for the coefficients of \i and v can be obtained for small a, performing a similar 
analysis for larger values of a is somewhat problematic. On that account, a numerical approach is followed instead. 
In particular, we propose to determine the appropriate values of fi and v such that the following error is minimized: 

oo 

£ = 6J2 {Qi(a, 5P) - exp (- e v {5f5Y)f (5) 

/3=0 

where 5 is small and £ — >• £ as <5 — >• 0. Such optimization problems can be solved using numerical techniques 
flU. It should be noted that the problem of minimizing the error term defined in §5§ is not a convex optimization 
problem. As such, numerical methods may not always converge to the global optimum. Nevertheless, we find that 
the observed optimum is often adequate, as illustrated in Section [TV] 

Since Qi(a,b) decays exponentially with b, we argue that we can ignore terms in the summation in ® 
corresponding to values of b larger than some 6 max in order to facilitate optimization. This is particularly justified 
by noting that we are interested in obtaining an accurate expression for Q\ that captures most of its mass. Thus, 
the upper limit on the summation in © can be replaced by /3 max = b maiX /5. As an example, we set 5 = 10~ 4 
and b max = 12, and solved the above optimization problem using a line-search algorithm for several values of a. 
Results are shown in Table U 

Using the values listed in Table H it is possible to derive an approximate expression for fi(a) and v{a) using 
polynomial regression J5J- For instance, assuming that p{a) is a polynomial of fourth order in a, the regression 
model for p,(a) can be expressed as 

n(aj) = Mo + + £2% + M3«| + /MOj + £j (6) 

for j = 1, . . . , N where ej is the error in the approximation, {£ti}f =1 are the estimation coefficients, and N is the 
number of observed instances (c.f. Table U). The above expression can be written in matrix form as /i = A/i + e, 
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Fig. 2. Plots of fi(a) and v(a) obtained through a numerical approach. The solid lines represent the optimized values, while the squares 
and circles show the approximate values from Q. 



where p, = [fiQ, ■ ■ ■ , /X4] T and /x is similarly denned, A is an iVx 5 matrix with kth column being 

for k = 1, . . . , 5, and e = [ei, • • • , ejv] ■ Using ordinary least squares estimation, the coefficients can be obtained 
using ft, = (A T A) 1 A T n. For m = 4, this approach yields 

fj,(a) = 2.174 - 0.592a + 0.593a 2 - 0.092a 3 + 0.005a 4 

u(a) = -0.840 + 0.327a - 0.740a 2 + 0.083a 3 - 0.004a 4 . (7) 

Fig. |2] depicts the comparison between the optimized values (from Table IJ) and approximated values of the two 
parameters given in f7]). It can be observed from the plots that the two sets of values are very close, indicating the 
suitability of the above two equations. Such approximations are convenient if fast computation of the parameters 
v{a) and n(a) are required. 



fr-l ... n*- 1 
1 ' ' N 



III. An Application of the Proposed Approximation 

As previously stated, the presented approximation can be particularly useful when the CCDF of the power of a 
Rician channel needs to be integrated over the second argument. The power of a Rician channel is noncentral-x 2 
distributed, whose CCDF is given by 

Fx (x) = Qi (V2K, ^/2u- 1 {K + l)x) (8) 

where K is the Rice factor and u> is a channel dependent parameter. Given that a = \llK in this case and, in 
general, 1 < K < 10 ifTOl . it follows that a < 5. On that account, the approximations of the parameters /x(a) and 
v(a) presented above can readily be used. 

To demonstrate the use of the proposed approximation of the Marcum Q-function, we consider the analysis of the 
full connection probability of a random ad hoc network, similar to the work presented in ifTTI . lfl~2l . Consider the 
connection probability of two nodes in a system, which we denote by H, given a minimum data rate requirement 
of Rq. By adopting an information theoretic definition of connectivity, we define 

H = P(log 2 (l + 1 \h\ 2 )>R ) (9) 

where \h\ 2 is the channel gain between the two nodes and 7 is the SNR which is dependent upon the distance 
between the two nodes and other parameters such as the path loss exponent and antenna gains. Under the assumption 
of a Rician fading channel with Rice factor K and a path loss exponent of two for illustration, we have H(r) = 
Qi (V%K, ra) where r is the distance between the two nodes and a is a function of the system parameters. To 
derive the probability that the network is fully connected, it is necessary to average H(r) over the configuration 
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e ub with approximate parameters 
Approximation in [7] 



with optimal parameters 
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Fig. 3. Comparison of the Marcum Q-function with the proposed approximation and that presented in |7|. For larger a, the approximation 
in (7) diverges from the actual curve for small b. Optimal parameters for the approximations are obtained by minimizing J5), while the 
approximate parameters are obtained from 0. 

space lfl2l . For a homogeneous system, this amounts to averaging H(r) over all distances between nodes. Such a 
calculation would involve an integral of the form (c.f., (19)-(21) in lfl2l for Rayleigh fading) 



where A = e u a^, r\ and r 2 are the minimum and maximum distances between nodes within the system, and 7(2;, y) 
is the lower incomplete gamma function. The complete analysis of the full connection probability is beyond the 
scope of this letter. What is important to note is that without the approximation derived in this paper, solving the 
integral stated above would be very challenging if not impossible. 



To demonstrate the accuracy of the proposed approximation, we compare the proposed method to an existing 
approximation in the literature [7] that generally yields small approximation error. The comparisons are shown 
in Fig. [3] For the approximation presented in Q, the value of k was set to 50 in equation (6) therein. It can 
be observed from the plot that, for small a, the proposed approximations are close to the Marcum Q-function. 
However, as a increases, divergence from the actual curve is seen for the approximation of (7 j. Nevertheless, the 
proposed approximations still adequately represent the mass of the Marcum Q-function over the range of b values; 
consequently, our approximation is robust with respect to changes in a. 

For large values of b, we note that existing bounds and approximations in the literature often provide a more 
accurate representation of Qi(a, b) compared to the proposed method. This can easily be observed graphically, but 
we omit the results here due to space constraints. For applications that require a closer approximation for large 
b, the expressions in H, Q and references therein would be more appropriate. However, if the integral of the 
Marcum Q-function over the domain of the second argument is sought, the approximation presented in this paper 
is more suitable. 

We next compare £ for our proposed approximation and the one given in Q. Results shown in Fig. [4] for different 
values of a demonstrate the accuracy of the proposed scheme. As mentioned in the previous section, the range 
of a values considered in the plot is the range that would typically be encountered in practice in communication 
system analysis. However, for the problem defined in ©, it is guaranteed that the solution would minimize the error 
term for any value of a. The reader should also bear in mind that the method presented in this paper minimizes 
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IV. Accuracy of the Approximation 
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Fig. 4. Comparison of the approximation error £ using the proposed approach and the one given in (7). The proposed approximation 
remains robust to changes in a.) 

a particular definition of the approximation error. Depending on the application, other objective functions may be 
defined. 

V. Conclusion 

In this paper, a simple approximation of the first order Marcum Q-function was presented that can be used 
in network connectivity analysis. For small input argument a, an analytical approach was presented for finding 
the approximation parameters, while for larger a, a numerical procedure based on an optimization problem was 
proposed. Equations for approximating these parameters were then presented. Simulation results demonstrated that 
the approximations led to an accurate representation of the Marcum Q-function, especially for small values of b. 
As b tends to infinity however, existing bounds of Q\(a, b) yield to a closer representation of the function. 
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